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Relation between the density-matrix theory and the pairing theory 

Mitsuru Tohyama and Satoshi Takahara 
Kyorin University School of Medicine, Mitaka, Tokyo 181-8611, Japan 

^j- ' The time-dependent density-matrix theory (TDDM) gives a correlated ground state as 

a stationary solution of the time-dependent equations for one-body and two-body density 
\ matrices. The small amplitude limit of TDDM (STDDM) is a version of extended RPA 

^s^j ■ theories which include the effects of ground state correlations. It is shown that the solutions 

of the Hartree-Fock Bogoliubov theory and the quasi-particle RPA satisfy the TDDM and 
STDDM equations, respectively, when only pairing-type correlations are taken into account 
in TDDM and STDDM. 

S: 

\Q ' §1. Introduction 

(N : 

', The quasi-particle random phase approximation (QRPA) based on the ground 

J> ' state in the Hartree-Fock Bogoliubov theory (HFB) has extensively been used as a 

standard microscopic theory to study nuclear structure problems 1 ^ see also Ref. 2 )~ 4 ) 
for recent applications). HFB and QRPA treat the pairing correlations in the frame- 
work of a mean-field theory: They are the stationary and small amplitude limits 
of the time-dependent HFB theory (TDHFB), respectively. On the other hand the 
time-dependent density-matrix theory (TDDM) 5 )' 6 ) treats two-body correlations, 
including the pairing correlations, as genuine two-body effects. TDDM gives a cor- 
related ground state as a stationary solution of the TDDM equations and the small 
amplitude limit of TDDM (STDDM) 7 ) becomes a version of extended RPA theo- 
ries which include the effects of ground-state correlations. This TDDM + STDDM 
scheme has recently been applied to nuclear structure problems. 8 ) Although the 
j> | TDDM + STDDM scheme deals with the pairing correlations, its relation to the 

HFB + QRPA scheme has not been investigated yet. The aim of this paper is to 
clarify a relation between TDDM and HFB for a correlated ground state and also a 
connection between STDDM and QRPA. We will show that when only pairing-type 
correlations are considered in TDDM and STDDM, the HFB and QRPA solutions 
satisfy the TDDM and STDDM equations, respectively. The paper is organized as 
follows: In sect. 2, the relation between TDHFB and TDDM is considered. The 
relation between the ground states in HFB and TDDM is discussed in sect. 3. The 
connection between QRPA and STDDM is given in sect. 4 and sect. 5 is devoted to a 
summary. 

§2. Relation between TDDM and TDHFB 

TDDM determines the time evolution of one-body and two-body density matri- 
ces p and p2 in a self-consistent manner. The equations of motion for p and p2 can be 
derived by truncating the well-known Born-Bogoliubov-Green-Kirkwood-Yvon hier- 
archy for reduced density matrices. 5 ) We expand p and C2, the correlated part of 
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p2, with a finite number of single-particle states ip a : 

p(ll',t) = ^n^(t)^(l)C'(l'), (2-1) 

C 2 (121 / 2',t)=p2-^(w) 

= CaWFWMVMWjVWpW, (2-2) 

where A(pp) is an antisymmetrized product of the one-body density matrices and the 
numbers denote space, spin and isospin coordinates. In the numerical applications of 
TDDM (for example Ref. 9 )) it is convenient to use a time-dependent single-particle 
basis which satisfies a TDHF-like equation. In the comparison between TDDM and 
TDHFB, however, a time-independent basis is appropriate. In such a basis the 
equations of motion of TDDM become the following two coupled equations: 6 ' 



ihh aal {t) = Y^aXUXa' - e\ a >n aX ) 
A 

+ [(«A 3 |i;|AiA2)C l AiA 2 a'A3 - CaA3AiA 2 (AiA 2 [i;ja / A3)], (2-3) 

A1A2A3 

ifoC a p a ip>(t) = "^^{e a \C\p a 'p' + ep\C a \ a ipi — e\ a iC a p\pi - epi\C a p a /\) 

A 

+ B a/ 3 a ipi + Papa'pi + Hapaip, (2-4) 

where v is a two-body interaction and e a p is the matrix element of the time-dependent 
mean-field Hamiltonian h defined as 

h(p)ip a (l) = --^—$JX) + / d2«(l,2)|>(22, t)ip a (l) -p(12,*)Va(2)]. (2-5) 

For simplicity we assume that v does not depend on the nuclear density. The term 
Bafia'P' ° n the right-hand side of Eq. (|2-4|) represents the Born terms (the first-order 
terms of v). The terms -P a /3a'/3' an d H a p a ipi in Ea. (|2-4I) contain C a p a ipi and repre- 
sent higher-order particle-particle (and hole-hole) and particle-hole type correlations, 
respectively. Thus full two-body correlations including those induced by the Pauli 
exclusion principle are taken into account in the equation of motion for C a j3 a ipi. 
The explicit expressions for B aj 3 a tp>, P a /3a>(3' and H a p a ipi are given in Ref. 6 ) and 
shown again in Appendix A for completeness. Equations ()2-3|) and (|2-4|) can also be 
obtained from 

ih^(${t)\a+a a \${t)) = ihh aa ,(t) = ($(t)\[a+a a ,H]\<l>(t)), (2-6) 

ih—($(t)\ : a+a+apda : \$(t)) = ihC a/3a > p (t) 

= ($(t)\[:a+a+apa a :,H]\<l>(t)), (2-7) 
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by expressing a three-body density matrix in terms of n aa >(t) and C a p a ip>(t). In 
the above equations, \$(t)) is the total wavefunction, H is the total hamiltonian 
consisting of the kinetic energy term and a two-body interaction, [ ] stands for the 
commutation relation, and : a^,a^,apa a : means 

: a^apapaa := a+,a^,apa a - A(ap,apn aa >(t) + a+,a a nppi{t)). (2-8) 

The equation for C a p a /pi(t) without the B and H terms has a relation with the 
solution in TDHFB. We write it explicitly: 

ihCafia' 0' (t) = y^^axCxpgip' + epxC a x a >p> — eXa'C a pxp> — ep'\C a p a 'x) + Pccpa'p' 
X 

= y^^axCxpgipi + epxCaXa'P' ~ ^Xa'C a pxp' ~ Cp< xC a f3a' x) 
X 

+ (A1A2IU \^3^4)[(S a X 1 SpX2 — 5aX i npx 2 — n aX 1 5pX 2 )Cx [i X i a'f3' 

A1A2 A3 A4 

— (^A 3 o'^A 4 /3' — £>A 3 a' n A 4 /3' ~ n A 3 a'£A4/3')Ca/3AiA 2 ]- (2-9) 

The TDHFB equations can be obtained from Ea. (|2-3j) and another equation for 
the pairing tensor n a p(t) = (<P(t)\apa a \<P(t)) 

ihk aP (t) = (<P(t)\[apa a ,H]\$(t)), (2-10) 

assuming that C a p a >p>{t) is a product of the pairing tensors, that is, C a p a /pi(t) « 
(<P(t)\a+,a+\<P(t))(<P(t)\apa a \<P(t)) = K* a rp,(t)K a p(t) and that ($(t)\a+,ap / apa a \$(t)) » 
A(n aa '(t)Kppi(t)). The wavefunction in TDHFB is different from that in TDDM. 
However, we use the same notation |#(i)) as long as it does not cause a confusion. 
Equation (|2-3|) becomes 

ihh aa <(t) = ^(eaxnxa' - ex a m a x) + ^J(A*A««'A ~ A* a , x K a \). (241) 

A A 

The pairing potential A aa /(t) is defined as 

A*a'(*) = 2 S (^'\v\X 1 X 2 )AKx 1 x 2 (t), (242) 

Ai A2 

where the subscript A indicates that the corresponding matrix is antisymmetrized. 
Equation (|2-10j) becomes 

ihk a p(t) = ^(eaxKxp + epx^ax) + A a p + ) (Ap\n a x - A aX np X ). (243) 

A A 
In TDHFB, Eqs. (|2-ll|) and (|2-13|) are solved under the constraint of the total number 



of particles introducing the Lagrange multiplier A and replacing e a p by e' a p = e a p — 
X8 a p. Using 7i and 1Z defined as 

«=(-!• -u) • < 2 - i4 » 
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where h' af3 = e' a/3 , we can express Eqs. (|2-ll|) and (|2-13|) and their complex conjugates 
in matrix form: 10 ) 

ifgR. = [H,H]. (2-16) 

The condition 

K 2 = K (2-17) 

holds due to the facts that the quasi-particle transformation is unitary and that 
|^(t)) in TDHFB is the vacuum with respect to the quasi particles. 1 )' 10 ) 

In the following we discuss the relation between the equations in TDDM and 
TDHFB. When the approximation C a p a tpi w K^ lj3 ,K a p is made, the first TDDM 
equation (Eq, (|2-3|) ) becomes Ea. (|2-ll|) and the second TDDM equations fEa. l|2-9|) ') 
is written as 

ih(k a/ 3(t)K* al/3l (t) + K a/3 (t)k* al/3/ (t)) 

= ^[(e Q A K A/3 + £/3A K aA)^a'/3' ~ ( e Aa' K A/3' + 6 A/3' K a> \) K af3] 
A 

+ [A a/3 - ^(A*a«/?a + A\pn a \)]K* a t p, 

A 

- [A* a , p , - + A* xf} ,n Xa ,)]K ap . (2-18) 

A 

Using Ea. (|2-13|) and its complex conjugate 

- ihk* a , p ,{t) = ^(e Aa 'K A( g/ + e A/3 'K*, A ) 

A 

+ A *«'P> + J^js'A'W " K^W), (2-19) 
A 

we evaluate the sum [Eq. (|2T3j) x k* — Eq. (I2T9|) x n a p\ . It is straightforward to 
show that the sum [Ea. lf^T^ x K* a ,m -Ea. (f2~T§)) x K a/3 ] is equivalent to Ea. (|2^TH|) . The 
replacement of e a p by e' a/3 in Eas. (|2T3j) and (|2T9() does not change the conclusion 
because there is a subtraction among the single-particle energies in Ea. (|2T8|) . Thus 
it is shown that the solution of TDHFB satisfies the TDDM equations fEas. (|2-3(l and 
(|2-9|) ) under the conditions that the correlation matrix C a p a ip> is constructed using 
n a p in TDHFB as C a p a i p «* ,q, K a p and that only the paring- type correlations are 
included in TDDM. 

§3. Relation between the HFB ground state and a stationary solution 

of TDDM 

The ground state |^o) m TDDM is a stationary solution of the TDDM equations, 
and the occupation matrix n^ a , = (^ol a a' a a|^o) an d the two-body correlation matrix 
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|^o ) are time independent. The equations for these 
quantities are obtained from Eas. (|23[l and (|2-4j) 

i A2C/A3 - C a\ 3 X 1 X 2 ( X l X 2\v\a A3}) 

A1A2A3 

= 0, (3-1) 

(e« + e/3 — e a i - ep^C^^^t + -^a/fo'/?' + Paj3a> P> + ^aPa'P' = °' ( 3 ' 2 ) 

where P°, and i?° are -B, P, and written in terms of n° and C°, respectively. 
Hereafter we use the eigenstates of the mean-field hamiltonian as the single-particle 
basis. We have shown in Ref. 8 ) that the solution of Eas. (|3T|) and (|3-2j) can be 
obtained using an iterative gradient method. Equation (|3-2|) without B° and H°, 
which has a relation with HFB, comes from Eq. ()2-9[) 

(e Q + ep- e a > - ep>)C a p a ,p, + P a p a/ p, = (e a + ep - e a > - ep>)C a p a ,p, 

+ Yl ( A 1 A 2M A 3A 4 )[(<W<W 2 -«5aAi^A 2 - n aM 5 pX2) C X^a'P' 
Ai A2A3A4 

- (^A 3 a'^A 4 /3' - ^A 3 a' re A4/3' ~ n W ^40' ) Ca/?Ai A3] = °- ( 3 " 3 ) 

The TDHFB equation Eq. Q2-ll|) becomes one of the HFB equations 

(e« - e a ,)n° aa , + £« a /$a - A% x K° aX ) = 0, (3-4) 

A 

where = {<Po\apa a \<Po) and A° a g is the pairing potential written in terms of «£g. 
This can also be obtained from Ea. (|3T|) assuming C£g a /g/ K ^p' K ap- Equation 
(|2- 13f) becomes the other HFB equation 

(e a + ep)n% + A% + £(4a<A - A° aX n o 0X ) = 0. (3-5) 

A 

We can express Eqs. (j3-4j) and (|3-5f) and their complex conjugates in matrix form: 1 )' 10 ) 

[H°,n ) = 0, (3-6) 

where TC° and 1Z° are TL and 1Z written in terms of n^ Q , and k^. The condition 
Ea. (|2T7|) becomes 

(ft ) 2 = n°. (3-7) 

In HFB, the quasi-particle basis which diagonalizes both Ti° and 1Z° is usually used. 

Although it may be obvious from the discussion in sect. 2 that the solution of 
HFB satisfies the TDDM equations (Eas. ()3-l|) and (|3-3|l ). we show it explicitly in 
the following. When the approximation C® g ,g, ~ i^wg/i\.0 * s m& de, the first TDDM 
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equation (Eq, (|3-l|) ) becomes Eq. (j3-4[) as mentioned above. The second TDDM equa- 
tion fEg. (|3-3jO becomes 

A 

- [<V - J2( A «'A' + A%,n%)]n% = 0.(3-8) 

A 

Using Ea. (|3-5|) and its complex conjugate 

A 

we evaluate [Ea. (|3~5)) x n°*,p, - Eq.(j33J x k°„], which becomes Ea. (|3~%|) . Thus it is 
shown that the solution of HFB satisfies the TDDM equations (Eqs. ([3T[l and (|3-3jl ) 

for stationary states. 

§4. Relation between QRPA and STDDM 

STDDM has been formulated by linearizing the TDDM equations. 7 ) The equa- 
tions of STDDM for the one-body transition amplitude x aa >(fi) and the two-body 
transition amplitude X a p a ip/(n) can be written in matrix form 11 ) 

(4-1) 

The matrices a, b, c, and d are given in Ref. 11 ) but shown in Appendix B for 
completeness. Ea. (j4-l|) can also be obtained from the following equations: 

(#o|[aJa a ,fll|#)=w M («Po|a>«|#> I (4-2) 
(<Z> |[: a+,apapa a :,J3]|#> = oj^o\ : a+a+apa a : \<P), (4-3) 

where |#) is the wavefunction for an excited state with excitation energy to^. Lin- 
earizing Eas. ()4-2j) and (|4-3(l with respect to x aa > = (^o|Oo' aQ l^) an< ^ X a p a >p> = (^o| : 
a a' a ~8' a P aa : l^)> we obtain Eq. (|4-l(l . To find a relation with QRPA, we explicitly 
write down the equations for x a p and X a p a > pi : 

-e a + e al )x aal = ^2 [((^2\v\X 3 a r }An aXl - (a\ 2 \v |A 3 Ai}An° ia ,)x A3A2 

Ai A2A3 

+ ^("Ai[i;jA2A3)A^A2A3Q'A 1 - ^(AlA 2 |u|a'A 3 )AA' Q A3A 1 A 2 ], (4-4) 

= Yl [(( aA il u l A 3 A 2>AC , ^ /3a , /3 , + (/3Ai|w|A 3 A2} J 4C , ° A3a , / 3, 

A1A2A3 

- (\ 3 \ 1 \v\a'\ 2 ) A C°p X3 p, - (\ 3 \ 1 \v\P'\ 2 )aC®p 01 , X3 )xx 2 \ 1 



a 










b 


2 




) =""( 
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+ ^(AiA 2 |t'|a / A 3 ) j4 C° /3AlA2 a; A 3 / 3/ + ~(AiA 2 MA 3 (3 ') AC% Xl x 2 x\ 3a > 

- ^(aX 3 \v\X 1 X 2 ) aC^^^x^ - ^(A3/3|^|AiA2)AC^ lA2Q , / 3,x a A 3 ] 

+ ^{[(a^|A 1 A 2 )^ 
A1A2 

- ^2({ aX 3\v\\i\2)An°p X3 + (A3/?|f|AiA 2 } J 4n° A3 )]X AlA2a / /3 / 
a 3 

- [{X 1 X 2 \v\a'0) A 

- ^((A^Ma'AaUnV + (A 1 A 2 |^|A3 ( 5 , )X 3Q 0]^a 1 a 2 }- (4-5) 

The equation for X a p a ipi given by Ea. (|4-l|) contains the terms which describe a gen- 
eral coupling to the one-body amplitudes as well as all kinds of two-body correlations. 
To make a comparison with QRPA, we consider only the pairing-type correlations 
in Ea. (|4-5j) : The terms corresponding to B and H in Ea, H2-4|) are neglected. 

The QRPA equations are usually formulated using the quasi-particle represen- 



tation but may also be obtained from the following equations 

(0 o \[a+,a a ,H]\0) = u tl ($o\a+,a a \$) = w lt x ac s, (4-6) 

(&o\[apa a ,H]\$) = w M (* |a/9Oa|*) = ^n^ap, (4-7) 

(0 Q \[ata+,H]\0) = Up{& o \a±a£,\0) = u^p. (4-8) 



Under the assumptions (0o\a+,apai3a a \0) w -4(n^,x QQ ' + n® aa ,Xf3p>) + K Q *,p,il) a p + 
K a/3<?W> (0o\ a t' a P' a P a o\®) ~ A(K°pp,x aa > + n° a a ,Tp/3f3> ), and (0 o \a^apa^a a \0) w 
A{K p lj} x aa i +n° a ,^ /3/ 3/), Eqs.(giSJ-(|^EIl become 

(a;,, - e Q + e a /)x Qa / = ^ [((aA 2 |u|AiA 3 )An AlQ , - (AiA 2 |t>|a'A 3 )An° Al )x A 3 A2 

Ai A2A3 

+ ^(aAi|v|A 2 A 3 )AK°T Al VA 2 A 3 + ^(AiA 2 |w|a'A 3 )AK° A 3^A 1 A 2 ] 
+ ^« a <Aa q '-<V/> q a), (4-9) 

A 

- e a - eptyap = ^(A^xpx ~ A° Xf3 x aX ) 

A 

+ \ [(«/?I«|AiA 2 )a + ^((/3A 3 | V |AiA 2 ) A < A3 + {\ 3 a\v\X 1 X 2 ) A np X3 )}^x 1 x 2 
A1A2 A3 

+ ^ ((X 1 a\v\X 2 X 3 ) A 4 3 p- (XMXiX^akIJx^, (4-10) 
A1A2A3 

(u^ + e a > + ep>)^p> a > = ^2(A° Xf3 ,x Xa i - A x * a ,x X p) 

A 
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- \ E [(AiA 2 |H/3'«% - y £((\ 1 \ 2 \vW\ 3 )An X3a , + (A 1 A 2 | V |A 3 a%<^)]0A 1 A 2 
A1A2 A3 

+ ^ ((A 1 A3|HA 2 a%4T A3 -(A 1 A3|^|A 2 /3%^ A3 ) 2; A 2 A 1 . (4-11) 
A1A2A3 

These equations correspond to the linearization of Ea. (j2-16Jl and are written in matrix 
form 3 )' 10 ) 

w^X = [H°,X] + [5H,K°], (4-12) 
where X and STL are defined as 

*-(;-$-)■ (443> 

W= |^ Ea 1 A 2 (« A iI u I« ,A 2) j 4Xa 2 A 1 iEa 1 A 2 (««>I A 1 A 2)a^A 1 A 2 \ 
V lEA 1 A 2 ( A l A 2b|a" / )A</'A 1 A 2 -EA 1 A 2 («' A ll t, l« A 2)A2;A 2 A 1 /' 

The linearization of Ea. H2-17j) becomes 

X = K°X + XTl . (4-15) 



The transformation of Eq. (|4T2|) using the quasi-particle basis and the conditions 
for 7Z° and X (Eas. l|3-7|) and fl|4- 15|t ) gives the QRPA equations for the amplitudes 
associated with the excitation operators (3^/3^ and fikftk,' 10 ' where p£ and /3k are 
the creation and annihilation operators of a quasi particle, respectively. 

In the following we make a comparison between QRPA and STDDM. Using 
X a pa'p' ~ Kap&P'a' + K a*/3'V'a/3 i n the equation for x aa ' (Eq. (|4-4|l ). we obtain the one 
of the QRPA equations Eq. (j4-9|h Similarly, the equation for X a p a ipi (Eq. ([4-5|) ) can 
be modified as 

(ay, - e a - ep + e a > + e >)(K O a p(f)p' a ' + K°*,p,i/) a p) 

= Yl [((a A l| u l A 3 A 2)AKA 3 /3 + (/ 3A l| u l A 3 A 2)AK°A 3 ) K °V 
AiA 2 A3 

- ((A3A1IV \a \ 2 ) AH%i3' + (A3Al|w|/3'A 2 )AKa*A 3 ) K a/3] a;; A 2 A 1 

+ E[^«'A X A/3' + A%,X Xa ,)K% - {Al x x px + A%X aX )K° a 1p,} 



+ - E« a ^a + A%nl x )]^ a , 

A 

- [Z^, - £(4$A»V + ^a>L')]^ 

A 

+ ^ 5^{[<«/3|i;|A 1 A 2 > A 



2 

AiA 2 



^((aA 3 |?;|AiA 2 )An° A3 + (A 3 /3|v |AiA 2 )An° A3 )]K°* /3 ,V'A 1 A 2 

As 
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+ [{X\Xi\v\a' 0) A 

- ^((A 1 A 2 |^|a / A 3 )An? A 3 /3 , + (A 1 A 2 | W |A 3 /3 , )A^ 3a ,)]^0A 1 A 2 }. (4-16) 

As 



Taking the su m [Eq . (l4~TTIl x K° a p - Eq. x + Eq. ff^TTT^ x K°Jp, + Eq. (fSZHJ) X ?M , 
we obtain Ea. (|TTB|) . Thus it is shown that the solutions of QRPA satisfy the STDDM 
equations (Eas. H4-4|) and (|4-5p ) when the two-body amplitude X a p a ipi is constructed 
from K Q a p in HFB, and ^ a/3 and </> Q/3 in QRPA as X a p a ipi rs K^^'a' + «°* / g/'0a/3- 

§5. Summary 

We investigated the relation between the time-dependent density-matrix the- 
ory (TDDM) and the time-dependent Hartree-Fock-Bogoliubov theory (TDHFB). It 
was shown that the two-body density matrix constructed from the pairing tensor in 
TDHFB satisfies the TDDM equation when only the pairing-type correlations are 
considered in the TDDM equation. As a natural consequence of special limits of 
the time-dependent theories, it was found that when only the pairing-type corre- 
lations are taken into account in TDDM and the small amplitude limit of TDDM 
(STDDM), the HFB and quasi-particle RPA (QRPA) solutions satisfy the TDDM 
and STDDM equations, respectively. Thus it is clarified that the HFB + QRPA 
scheme gives a subset of the solutions of the TDDM+STDDM scheme where the 
two-body correlation matrix and transition amplitudes are of separable form. 

Appendix A 
A 

The terms B, P, and H in Eq. (|2-4|) are shown below: 



B a /3a'f3> = (AlA 2 |w|A 3 A4) J 4[(5aAi - «qAi)(^/3A 2 _ n /3A 2 ) n A3a'«A 4 /3' 

A1A2A3A4 

- n aXl np\ 2 (d\ 3 a> ~ n\ 3a > ) (8\ A pi - nx 4 p< )] , (A- 1) 

Paf3a'f3' = (Al A 2 | V | A3 A4} [(5 a \ 1 5p\ 2 — 5 a x 1 U/3X 2 — n Q Ai <5/3A 2 )Cx s X 4 a'/3' 

A1A2A3A4 

— {°~X 3 a'°~X 4 l3' — 0~\ 3 a>nx 4 l3> ~ n X 3 a'$X 4 l3> )Ca/3X 1 A 2 ] j (A-2) 

H a /3a'j3' = {Xl\2\v\\3X4)A[8a\ 1 (nx 3 a'Cx 4 f3\ 2 l3' — nx 3 f3>Cx 4 f3X 2 a') 

A1A2A3A4 

+ 0~l3\ 2 { n \ 4 l3'C'X [i a\ 1 a' ~ n X 4 a' Cx^aX!^ ) 

— <Va 3 ( re aAi Ca 4 /3A 2 /3' ~ ^/3AiCA 4 aA 2 /3') 

- VA4( n /3A 2 CA3aAia' ~ n a x 2 Cx 3 /3X ia ')]- (A-3) 
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Appendix B 
B 

The matrices in Ea. (|4-lj) are shown below. The terms in b and d which are 
considered in Ea. (|4-5|) are underlined. 

a(aa' : AA') = (e Q - e a ,)5 a x6 a > X ' 

- ^2((f3\'\v\a'\) A n% - (a\'\v\(3X) A n$ a/ ), (B-l) 

P 

b (aia^Q^a^ : AA') 

= -<W{^((<W - n ° a 2/3) n ya' 1 n °5 a ' 2 + ™° 2 /3(<Vi ~ ^X'W, ~ ™L 2 ) <A'/?| v|-y«5> A 

Pi 

PlS 

/3 7 

+ ^A'-Q^fOW, - "L^^^^T + n L,(<W - "-° 1/ 3)(^27 - n ° 27 )(/ ? 7b|A5)A] 
PlS 

Pi 

- *q' 2 v{5ZK 5 K ~ n L;) n °i/3 n a 27 + - "ai/OOW - "-° 27 )(/37b|A5)A] 

/3 7 6 

+ Et Wj^llg^r + <"i^IHA7)aC° 2Vi/3 - («2/?b|A 7 )AC° i7a , i/3 ]} 

/3 7 

+ ^ [(aiA>|/?A) A C° a2a , a , - (a 2 A>|/?A) A C° aia , a , 

- Wj^Aj^^ + {pX'lvla'^AC^J , (B-2) 

c(aa' : AiA^A^) = (aA 2 |u|AiA 2 )(5 Q ,/y i - (A^A^u |a'A 2 )<5 QAl , (B-3) 

d(aia 2 a' l a l 2 ■ AiA 2 A' 1 A / 2 ) = (e ai + e a2 - e a / - e^^Ai^A^^^A!, 

+ ^A'/c^Af, Yl 0Wfgg7 ~ ^2 7 n °!/3 ~ <W^° 27 )(/?7MAiA 2 ) 

Pi 

- 5 aiXl 5 a2 x 2 ^2 (<5 Q ' i/3 <5 Q27 - S a ' 27 n° a , i/3 - S^n^iX'^v]^) 

Pi 
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J2((^'i\v\(3\i)An° pali 

P 



P 



(fi^2\v\ot 2 \i)An ail3 



+ ^aiAi^A'j 



(PX 2 \v\a' 2 X 2 ) A n o a20 



^(<a 2 A>|/3A 2 ) A n° , 

P 



(f3\' 1 \v\a' 1 X 2 ) A n° ol2 ^ 



(B-4) 



References 



1) P. Ring and P. Schuck, The Nuclear Many-Body Problem, (Springer- Verlag, Berlin, 1980). 

2) M. Matsuo, Prog. Theor. Phys. Suppl. 146 (2002), 110. 

3) E. Khan, N. Sandulescu, M. Grasso and N. Van Giai, Phys. Rev. C66 (2002) 024309. 

4) N. Yamagami and N. Van Giai, Phys. Rev. C69 (2004), 034301. 

5) S. J. Wang and W. Cassing, Ann. Phys. 159 (1985), 328; W. Cassing and S. J. Wang, Z. 
Phys. 328 (1987), 423. 

6) M. Gong and M. Tohyama, Z. Phys. A335 (1990), 153. 

7) M. Tohyama and M. Gong, Z. Phys. A332 (1989), 269. 

8) M. Tohyama, S. Takahara and P. Schuck, nucl-th/0311016 Euro. Phys. J. A in press. 

9) M. Tohyama and A. S. Umar, Phys. Lett. B549 (2002), 72. 

10) J. G. Valatin, Phys. Rev. 122 (1961), 1012. 

11) M. Tohyama and P. Schuck, Euro. Phys. J. A 19 (2004), 215. 



